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ABSTRACT 

Context. The Sun has recently been predicted to be an extended source of gamma-ray emission, produced by inverse-Compton (IC) 
scattering of cosmic-ray (CR) electrons on the solar radiation field. The emission was predicted to be extended and a confusing 
foreground for the diffuse extragalactic background even at large angular distances from the Sun. The solar disk is also expected to 
be a steady gamma-ray source. While these emissions are expected to be readily detectable in the future by GLAST, the situation for 
available EGRET data is more challenging. 

Aims. The theory of gamma-ray emission from IC scattering on the solar radiation field by Galactic CR electrons is given in detail. 
This is used as the basis for detection and model verification using EGRET data. 

Methods. We present a detailed study of the solar emission using the EGRET database, accounting for the effect of the emission 
from 3C 279, the moon, and other sources, which interfere with the solar emission. The analysis was performed for 2 energy ranges, 
above 300 MeV and for 100-300 MeV, as well as for the combination to improve the detection statistics. The technique was tested on 
the moon signal, with our results consistent with previous work. 

Results. Analyzing the EGRET database, we find evidence of emission from the solar disk and its halo. The observations are compared 
with our model for the extended emission. The spectrum of the solar disk emission and the spectrum of the extended emission have 
been obtained. The spectrum of the moon is also given. 

Conclusions. The observed intensity distribution and the flux are consistent with the predicted model of IC gamma-rays from the 
halo around the Sun. 

Key words. Methods: data analysis, Moon, Sun: X-rays, gamma rays, Gamma rays: observations, Gamma rays: theory 



(N 1. Introduction 



Solar quiet gamma-ray astronomy started playing a sign ificant 
role in the early 90's thanks to the EGRET mission. iHudson I 
dl989l) pointed out the importance of gamma-ray emission from 
the quiet Sun as an interesti ng possibility for h ighly sensitive 
instruments such as EGRET. Seckel et al.l (1 199 lb estimated that 
the flux of gamma rays produced by cosmic-ray interactions on 
the solar surface would be detectable by EGRET. This emission 
is the result of particle cascades initiat ed in the solar atmosp here 
by Galactic cosmic rays. Meanwhile Thompson et al.l d 19971) de- 
tected in the EGRET data the gamma rays produced by cosmic 
ray interactions with the lunar surface. Although they expected 
similar interactions on the Sun, the EGRET data they analyzed 
did not show any excess of gamma rays consistent with the po- 
sition of the Sun. They obtaine d an upper limit of the flux above 
100 MeV of 2 x 10" 7 cirr 2 s -1 . iFairbairn etafl d2007l) . analyzing 
the EGRET data of the solar occultation of 3C279 in 1 99 1 , found 
a photon flux of about 6 x 10~ 7 crrr 2 s~' above 100 MeV in the 
direction of the occulted source, which they used to put limits on 
axion models; however they did not consider the extended solar 



emission. 



Inverse Compton scattering of cosmic-ray electrons on the 
solar photon halo around the Sun has bee n estimated to be 
an important source of gamma-ray emission dOrlando & Strong] 



|2007at|Moskalenko et al. 2007)fl In our previous work, we pre- 
dicted this to be comparable to the diffuse background even at 
large angular distance from the Sun. The formalism has now 
been improved and more accurate estimates have been obtained. 
The anisotropic scattering formulation and the solar modulation 
have been implemented, in order to give a more precise model 
for the EGRET data. In this work we explain in detail our model 
of the extended emission produced via inverse Compton scatter- 
ing of cosmic -ray electrons on the solar radiation field. 

A first report of o ur detection with EGRET data was given in 
lOrlando et al.l d2007l) . Here we describe the full analysis includ- 
ing the spectrum of disk and extended components. Our result is 
very promising for the forthcoming GLAST mission which cer- 
tainly will be able to detect the emission even at larger angular 
separation from the Sun. Moreover, the extended emission from 
the Sun has to be taken into account since it can be strong enough 
to be a confusing background for Galactic and extragalactic dif- 
fuse emission studies. 

2. Theoretical model 

Inverse-Compton scattering of solar optical photons by ~ GeV 
cosmic -ray electrons produces y-radiation with energies of 100 
M eV and above. Improving o n the approximate estimates given 
in lOrlando & Strond (|2007a), in this paper the anisotropic for- 
mulation of the Klein-Nishina cross section has been introduced, 



Send offprint requests to: elena.orlando@mpe.mpg.de 



1 The same p r ocess fro m nearby stars has been studied in 
lOrlando & Strong i2007ah and lOrlando & Strong! i2007bh 
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Fig. 1. Definition of variables for eq.© describing the geometry 
of the solar inverse Compton emission. 



taking into account the radial distribution of the photons propa- 
gating outward from the Sun. Our elegant analytical formulation 
has been replaced by a numerical solution depending on the an- 
gle between the momenta of the cos mic ray electrons and the 
incoming photons. As pointed out in lOrlando & Stron g (2007a) 
this will affect the results at about the 10% level. Moreover, 
while in our previous work we used the modulated cosmic-ray 
electron spectrum as observed at earth, here a formulation of the 
solar modulation as a function of the distance from the Sun has 
been taking into account. 

The gamma-ray intensity spectrum is: 



KEy) 



4tt / 



e(Ey)dx 



(1) 



where the emissivity e is given by: 



e(E 7 ) = dE e x 



x J cr K _ N (y, E ph , E y ) n ph (E ph , r) c N(E e , r) dE ph 



(2) 



N(E e , r) is the electron density, with E e electron energy, 
n p h(E p i,,r) the solar photon density as function of the distance 
from the Sun r and the solar photon field E p h, cr is the Klein- 
Nishina cross section and y = E e lm e . The cosmic-ray electrons 
have been ass umed isotropically distrib uted everywhere in the 
heliosphere. In Orl ando & Strong! d2007al) . an elegant analytical 
formulation in eq © was obtained for the isotropic case, assum- 
ing a simple inverse square law for the photon density for all 
distances from the Sun and a cosmic ray electron spectrum con- 
sidered constant for all distances and determined from experi- 
mental data (see Fig.©) for 1 AU. In order to obtain the inverse 
Compton radiation over a line of sight at an angle a from the 
Sun, the photon field variation over the line of sight has to be 
known. As shown in Fig.©, for a given point x on the line of 
sight, the surface photon density is proportional to 1 /r 2 , where 
r 2 = x 2 + d 2 - 2 x d ■ cos a, with d distance from the Sun. 
Integrating the photon density over the line of sight from x — 
to x = oo one obtains: 



npi,(Eph) R 



dx 



2 (d ' cos a ~ x \ 

— —n P h\Eph) R arctan — /(a ■ sin a) — 

\ d- sin a ) 

2 I n/2 + arctan(cot a) 



n p h{E p h) R 



d ■ sin a 



(3) 



which depends only on the angle a, where R is the solar radius 
and HpkiEph) is the solar photon density. For large distance from 
the source n / ,/,(£ , p /,)=l/4 nBB(E p i,), where nBB(E p h) is the black- 
body density of the Sun. Integrating over solid angle and using 
eq.© and eq.© with the Klein-Nishina cross section, the total 
photon flux produced by inverse Compton scattering within an 
angle a becomes: 

i r 2n r r 

I(Ey) — — I d(p I sin a da | dE,,i, 
4?r Jo Jo J 

x J (T K - N {y, E ph , Ey) c N(E e ) dE e n ph {E ph ) R 2 x 

f 

Jo 



x 2 + d 2 — 2 x d ■ cos a 



16 d 



na + i^) 2 



■ arctan (cot 



a)) 



x J dE ph J cr KN c N(E e ) n BB (E p h) dE e (4) 

which for small a is proportional to a/d and the intensity / 
(per solid angle) is proportional to l/(ad). For the case of the 
anisotropic formulation of the Klein-Nishina cross section and 
for the electron modulation along the line of sig ht, we have to 
use nu merical computations, already adopted in Orland oet al.l 
(2007). 



2.1. Solar photon field 
2.1.1. Basic relations 

The Sun is treated as black body where the energy density 
is characterized by the effective temperature on the surface 
(T=5777 K) following the Stephan-Boltzmann equation. For the 
photon density close to the Sun, the simple inverse square law is 
inappropriate. In this work the emission from an extended source 
has been evaluated. The distribution of photon density from the 
Sun, as extended source, is given by integrating over the solid 
angle with the variables shown in Fig.© 




Fig. 2. Variables involved in eq.© for the calculation of the 
photon density around the Sun, where R is the solar radius and r 
is shown also in Fig. © 
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= 0.5 n BB (E ph ) [l - Vl " (R/r) 2 ] 



(5) 



with cos p=(r cos <p - R)/s and cos 4>max - R/r. For large 
distances from the Sun it reduces to the inverse square law. 

2.1 .2. Deviations of solar spectrum from black body 

We checked the effect of deviations of the solar spectrum from 
a black body. We used the black-body approximation with the 
AMO solar spectrum Q, referring to 2007, period of solar mini- 
mum, extrapolated back to the Sun. Since the estimated inverse 
Compton emission from the Sun is not affected significantly 
(<1%) considering the actual spectrum, we decided to keep the 
black body approximation for simplicity. We also verify that the 
EUV and XUV range of the solar spectrum do not affect the 
inverse-Compton emission. The irradiance at the shortest wave- 
lengths varies significantly with solar conditions and the range of 
variability spans a few percent at the longest wavelengths, from 
40% to more than 500% across the EUV and XUV , from solar 
minimum to solar maximum. Moreover, during solar flare, XUV 
increases, reaching factors of eight or more, and EU V increases 
on the order of 10 -40% dEparvier & Woods! 12003). However, 
lOncica et alJ d2002) found that the X-ray flares cannot contribute 
to the total solar irradiance fluctuations and even the most ener- 
getic X-ray flares cannot account for more than 1 mW/ra 2 . We 
found that the estimated inverse-Compton emission from the Sun 
is not affected significantly by the XUV/EUV flux and its vari- 
ability. 

2.2. Comparison of isotropic/anisotropic formulations 

In lOrlando & Strong! d2007al) the isot ropic formulation for 
the Kl ein-Nishina cross section given in iMoskalenko & Strong! 
(2000) was used in the analytical formulation for estimating the 
IC gamma-ray emission . In the present analysis we used the fol- 
lowing anisotropic Klein-Nishina cross section i n a more conve- 
nient form than in Moska lenko & Stron g (2000): 



0~K-N(y,E p h,E y ) 



E nhE 2 



F' 



1-v E'(l-v) 



+ (l-v) + 



1-v 



where v = E y /E e and y=E e /m e , with m e electron mass and 
E'ph = jE p h(\ + cos if) 



(6) 



(7) 



where rj is the scattering angle for the relativistic case, with r\ = 
for a head-on collision. In Fig.® is shown the contribution to 
the flux as a function of distance along the line-of-sight for the 
isotropic and anisotropic cross-sections. 

The integrated inverse-Compton emission was compared for 
isotropic and anisotropic cross sections; the maximum difference 
of the intensity for the two Klein-Nishina cross sections is of 
about 15% around 300 MeV and for 0.5° angular distances from 
the Sun. 

2.3. Electron spectrum 

The solar modulation of cosmic ray electrons has been 
widely studied, but is still not completely understood. 




Fig. 3. Relative contribution (arbitrary units) of the isotropic 
(gray lines) and anisotropic (black lines) cross sections, multi- 
plied by the photon density along the line of sight for E c =10 4 
MeV and E y =100 MeV. x and d are defined in Fig 1. The lines 
represent different angles from the Sun, from 1° to 6° (right to 
left). 



taken from http://www.spacewx.com/ 



Muller-M ellin et all d 1977b analyzed the data of Helios- 1 from 
1974 to 1975, a period that corresponds to the approach to so- 
lar minimum, and reported the measurements of proton and 
helium cosmic ray intensities for the 20-50 MeV/n energy in- 
terval between 0.3 and 1 AU. The proton gradients they ob- 
tained were small and generally consistent with zero, while 
the helium gra dients were positiv e but small. On the basis of 
Helios-1 data, dKunow et al.l ll975) stated: 'The radial variation 
seems to be v ery small for all pa rticles and all energy ranges'. 
The study of (Seckel et al. 1991) was based on this hypothe- 
sis. Since that time, the modulation theories have been devel- 
oped including effects lik e drifts and current sheet tilt variations 
(Joki pii & Thomas| [l981) and short time variations. In this work 
we care only about long period modulations for solar maximum 
and minimum and we assume that the longitude gradient is negli- 
gible. (McDonaldietal](|2003), analyzing cycle 20 and 22 of solar 
minimum between 15 and 72 AU, found that most of the modu- 
lation in this period occurs near the termination shock (assumed 
to be at 100 AU). On the other hand, from solar minimum to so- 
lar maximum the modul ation increases mainly inside the termi- 
nation shock. Recently Mo rales-Olivaresi & Caballero-Lopezl 
(120071) . investigating the radial intensity gradients of cosmic rays 
from 1 AU to the distant heliosphere and interpreting the data 
from IMP8, Voyagers 1 and 2, Pioneer 10 and BESS for Cycles 
21, 22 and 23, foun d different radial gradien ts in the inner helio- 
sphere compared to [McDonald et all d2003). In this region they 
obtained an average radial gradient of 3% /AU for 175 MeV 
H and =* 2.2% /AU for 265MeV/n He, that , at 1AU, give an in- 
tensit y smaller than [McDonald et al.l (120031) for H and higher for 
He. I Gieseler et all (120071) analyzed the data from Ulysses from 
1997 to 2006 at 5AU. They found a radial gradient of 4.5%/AU 
for a particles with energies from 125 to 200 MeV/n, which is 
consistent with previous measurements. 

In this work, our model of the modulation of cosm ic rays 
uses the formulation given in iGleeson & Axfordl |l968), using 
the studies regarding the radial distribution of cosmic rays in 
the heliosphere at solar minimum and maximum. The cosmic 
ray electron spectrum is given by the well known force field 
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approximation of cosmic ray nuclei used to obtain the mod- 
ulated differen tial intensity J(r,E) at en ergy E and distance r 
from the Sun (jGlee son & Axford| [l968l) . As demonstrated by 
Caballero-Lopez & Moraal ( 2004) it is a good approximation for 
Galactic cosmic rays in the inner heliosphere. The differential 
intensity of the modulated spectrum is given by: 

E(E + 2E ) 

J(r,E) = J(oo, E + <t>(r)) (8) 

(E + 0>(r) + 2Eq)(E + <D(r)) v ; 

where E is kinetic energy in MeV and Eq is the electron 
rest mass. We use the local interstellar electron spectrum 
J(oo, E + <£(>)). O(r) = (Ze/A)4> where <p is the modulation 
potential, Z the charge number and A the mass number. In order 
to compute the modulation potential i n the inner heliosphere, we 
used the parameterization found by Fuiii & McDonald! (2005) 
for Cycle 21 from 100 AU to 1 AU, neglecting the time depen- 
dence and normalizing at 1 AU for 5 00 MV (solar minimum) 
and 1000 MV (solar maximum). As iMoskalenko et alj d2007l) 
we find: 

(5(r) = <D (r- ai - r b 0A )/(l - r° b A ) (9) 

where (Do is the modulation potential at 1 AU, of 500 and 1000 
MV for solar minimum and maximum respectively; r/, = 100AU 
and r is the distance from the Sun in AU. Figure © shows the 
local interstellar electron spectrum ( Strong et al. 2004a) and the 
modulated spectrum at 1 AU, compared with the data. The elec- 
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Fig. 4. Loc a l inters tellar electron spectrum (blue line) as in 
IStrong et all d2004al) and the modulated at 1 AU of 500 MV 
(re d line) and 1000 MV (black line) compared with the data. 
See lStrong et al.1 (l2004al) for data reference. 

tron intensity as a function of distance from the Sun is given in 
Fig © for different energies and for modulation potential 500 
MV. 

The major contribution to the energy range 100-500 MeV 
we are interested in comes from electrons with energy between 
1 and 10 GeV, as shown in Fig. ©, obtained for the local inter- 
stellar electron spectrum, with no modulation. To compute the 
inverse-Compton extended emission from the heliosphere and 
compare it with the EGRET data, details of the solar atmosphere, 
the non-isotropic solar wind and the asymmetries of the mag- 
netic field have been neglected, considering the limited sensitiv- 
ity of EGRET. Since the biggest uncertaintes in the gamma-ray 
emission come from the cosmic -ray electron spectrum close to 




distance (AU) 

Fig. 5. Electron intensity as a function of distance from the Sun 
in AU, where corresponds to the Sun, for different energies and 
for modulation potential <t> = 500MV. 




Energy (MeV) 

Fig. 6. Contribution to the IC emission for different electron en- 
ergy ranges, left to right: 100 - 10 3 MeV, (1 - 3) x 10 3 MeV, 
(3 - 5) x 10 3 MeV, (0.5 - 1) x 10 4 MeV,(l - 5) x 10 4 MeV, 
(0.5 - 1) x 10 5 MeV, 1 x 10 5 - 5 x 10 7 MeV. 



the Sun, in our model we considered two possible configurations 
of the solar modulation within 1 AU. The "naive" approxima- 
tion is to assume that the cosmic -ray flux towards the Sun equals 
the observed flux at Earth, since there is evidence that modula- 
tion by the solar wind does not significantly alter the spectrum 
once cosmic rays have penetrated as far as Earth. Moreover, high 
energy electrons are less sensitive to the modulation. This ap- 
proximation gives an upper limit of the modelled flux. The other 
approach is to assume that the electron spectrum varies due to 
solar wind effects within 1 AU. With this "nominal" approxima- 
tion we assume that the formulation for solar modulation from 
100 AU to 1 AU can be extrapolated also below 1 AU, using eq. 
©. This gives an approximate lower limit in our model. 

2.4. Calculated extended solar emission 

Figure © shows the spectrum of the emission for two different 
angular distances from the Sun, without modulation and for two 
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levels of solar modulation ((D=500, 1000 MV, respectively for 
solar minimum and solar maximum) and for the cases of upper 
and lower limits described above. 




* ~* j„3 _,„5 _„6 

10 10 10 10 10 10 

Energy (MeV) 



Fig. 7. Spectrum of the emission for (top to bottom) 0.5° and 
5° angular distances from the Sun and for different conditions 
of solar modulation. Solid lines: no modulation, pink lines: 
<£o=500MV, blue lines :<E>o=1000MV, dashed lines: naive model, 
dotted lines: nominal model. EG B is the extragalactic back- 
ground as in lStrong etail (l2004bl) 

The angular profile of the emission is shown in Fig.® above 
100 MeV without modulation, for two levels of solar modula- 
tion (O=500, 1000 MV) and for the cases naive and nominal. 
The emission is extended and is important compared to the ex- 
tragalactic background (around 10 _5 cm -2 s -1 sr -1 ) even at very 
large angles from the Sun. Even around 10° it is still about 
10% of the extragalactic background and with the senstivity of 
GLAST should be included in the whole sky diffuse emission. 
An example of the intensity of the inverse-Compton emission 



E>100 MeV 




Angular distance from the Sun (degrees) 



Fig. 8. Angular profile of the emission as a function of the an- 
gular distance from the Sun above 100 MeV. Blue line: no mod- 
ulation, pink lines: <l>o=500MV, green lines: <E>o=1000MV, solid 
lines: naive model, dashed lines: nominal model. 

predicted in a region of 10° from the Sun is shown in Fig.® 



for 300-500 MeV. While the inverse Compton emission is ex- 




Fig. 9. Inverse-Compton emission modelled for a region of 10° 
from the Sun for 300-500 MeV and for the naive model of 1000 
MV modulation. Intensity is given in ran -2 s -1 sr _1 . In the figure, 
the minimum value of the intensity is 1.2 x 10~ 7 cm^s^sr -1 . 

pected to be readily detectable in future by GLAST, the situation 
for available EGRET data is more challenging. In the following 
section we present our study with the EGRET database. 

3. EGRET data preparation and selection 

We analyzed the EGRET dat a using the code developed for the 
moving target such as Earth (Petrv 2005). The software permits 
us to produce images centred on a moving source and traces of 
other sources around with respect to the centred source. To per- 
form the analysis for the Sun, we analyzed the data in a Sun- 
centred system. In order to have a better sensitivity and hence 
detection of our emission component, the diffuse background 
was reduced by excluding the Galactic plane within 15° lati- 
tude. Otherwise all available exposure from October 1991 to 
June 1995 was used (viewing periods 1, 11, 12, 19, 26, 209, 
221, 320, 325, 410, 420), excluding the big solar flare in June 
1991. When the Sun passed by other gamma-ray sources (moon, 
3C 279 and several quasars), these sources were included in the 
analysis. In order to obtain the flux of the moon to be fitted, the 
same method performed for the Sun h as been made with all data 
moon centred. Details will be given in lPetrv I (f2008). 

4. EGRET analysis 

We fitted the EGRET data in the Sun-centred system using 
a multi-parameter likelihood fitting technique (see Appendix). 
The region used for fitting is a circle of radius 10° around 
the Sun. The total number of counts predicted includes inverse 
Compton and solar disk flux, the moon, QSO 3C279, other back- 
ground sources and the background. We left as free parameters 
the solar extended inverse-Compton flux from the model, the 
solar disk flux (treated as point source since the solar disk is 
not resolvable by EGRET), a uniform background, and the flux 
of 3C279 - the dominant background point source. The moon 
flux was determined from moon-centred fits and the 3EG source 
fluxes were fixed at their catalogue values. All components were 
convolved with the energy-dependent EGRET PSF 
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Table 1. Parameters of the background sources used for the anal- 
ysis. 



Source 


Flux (> 100 MeV£ 


spectral index 


J0204+1458 


2.36 


2.23 


J0215+123 


1.80 


2.03 


J 1230-0204 


1.13 


2.85 


J1235+0233 


1.24 


2.39 


J 1246-0651 


1.29 


2.73 


J1310-0517 


1.04 


2.34 


J 1409-0745 


2.74 


2.29 


J2321-0328 


3.82 


_ b 



a Units 10~ 7 cirrV 

* The choice of the spectral index we used is explained in the text 



Table 2. Fitting results for the moon 



Energy (MeV) Moon flux (cm 2 s l ) Background (cm 2 s l sr ') 

TWO (5.55 ±0.65) xl0~' 3.47 xlO" 3 

> 300 (5.76 ± 1.66) xl0~ 8 1.01 xl0~ 5 

100- 300 (4.98 ±0.57) xl0~ 7 2.13 xl0~ 5 



4.4. Diffuse background 

Since the count maps used for our analysis include only data 
above |b| >15° latitude, and the Galactic plane is excluded, the 
Galactic emission can be well approximated by an isotropic 
background. We left its intensity as a free parameter in the fit- 
ting analysis. 



4.1. Model of extended solar emission 

The solar extended inverse-Compton emission was modelled 
within 10° radius from the Sun and with a pixel size of 0.5°. 
The model was convolved with the energy-dependent EGRET 
PSF and implemented in the fit, with a scale factor as free pa- 
rameter. We took into account the different approximations to 
the modulation described above, since the data we used cover a 
period from October 1991 when there was a solar maximum to 
June 1995 close to the solar minimum. Since the emission we 
are looking for is very close to the EGRET sensitivity limit, we 
had to take as much exposure as possible and it was not possi- 
ble to perform an analysis splitting the data according to solar 
modulation. 



4.5. Moon 

The flux of the moon was determined from moon-centred data, 
fitting the flux and an isotropic background. Since the moon 
moves quickly across the sky, all other sources moving by are 
taken into account as a constant component included in the 
isotropic background. The fit was performed for 2 energy ranges: 
100-300 MeV and >300 MeV. The value of flux of the moon 
for the different energy ranges was obtained by maximum likeli- 
hood. The maximum likelihood ratio statistic used for this anal- 
ysis is described in the appendix. The values of the best fit fluxes 
and lcr errors are given in Table ||2). Previous studies of the 
EGRET data (Thompson et al.) gave a flux of (5.4 ± 0.7) x 10~ 7 
cnr 2 s~' above 100 MeV in agreement with our analysis. The 



4.2. 3C279 

Since 3C 279 is the brightest source that passes close to the 
Sun, we decided to leave its flux as free parameter in the fit. We 
fixed the spectral index at 1 .96, the value given in the EGRET 
catalogue. The flux obtained was then compared to that in the 
EGRET catalogue. 

4.3. Other background sources 

Many sources passing within 10° from the Sun were included 
in the analysis. We decide d to implem ent their fluxes as fixed 
parameters. With the code (Petry 2005) it was possible to know 
the exact viewing period of any source a nd hence to determin e 
its flux from the 3rd EGRET catalogue dHartman et al.lll999h . 
Hence, we took the flux value in the EGRET catalogue corre- 
sponding to the period we used, when this was listed. Since the 
catalogue does not contain the flux values corresponding to all 
the periods in which the sources were close to the Sun, in case 
of incertainty we decided to use the first entry in the catalogue, 
which is the one from which the source position was derived. 
In almost all cases, this is the detection with the highest statis- 
tical significance, which in most cases corresponds to the max- 
imum value of the flux. We also took the spectral index from 
the EGRET catalogue. The values adopted for the sources are 
listed in Table (fl}. Since the spectral index of J232 1-0328 is 
not known, the fit was performed without this source above 300 
MeV, since in the literature there is no evidence of its detection 
above this energy. In order to verify that this choice does not af- 
fect significantly the fit result, we tested also a typical spectral 
index of 2. Both fitting results will be reported in the following 
results section. The traces of the background quasars were added 
together in the same data file, rescaled according to their fluxes. 















p 



10 3 

Energy (MeV) 



Fig. 10. Spectrum of the moon obtained by the fitting anal- 
ysis of the EGRET data between October '91 to June '95 
(dashed region). Thickness of the line shows the error bars. 
Black bars are obtained by Thompso n et alJ d 19971) for differ- 
ent solar conditions, while lines are the theoretical model of 
Moskalen ko & Porter! (120071) . for solar maximum (red line) and 
solar minimum (green line). 

spectrum of the moon shown in Fig.dTOb has been obtained im- 
posing a constant spectral index y, extrapolated up to 2 GeV 

The lunar spectral index was determined as 3.05*Q2g. Errors 
are calculated from the uncertainties on the integrated fluxes. 
Since the EGRET data we used to analyze the emission from the 
Sun extend from solar maximum to minimum, in this analysis 
we estimated an average lunar flux between the two periods of 
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solar modulation. Comparin g our spectrum with the model of 
Mos kalenko & Ported ([2007) for different solar conditions, we 
find a good agreement at 100 MeV, while at higher energies the 
model is about a factor of tw o higher. More extensive analysis 
will be given in lPetrvl d2008l) . 

5. Solar analysis results 

The analysis was performed for 2 energy ranges : 100-300 MeV 
and above 300 MeV as well as the combination to improve the 
detection statistics. This choice was determined by the limita- 
tions of the EGRET data. FigQj] shows the smoothed count 
maps, centred on the Sun, for those energy ranges and for a 
region of 20° side. The fitting region has a radius of 10° cen- 
tred on the Sun. Since the interesting parameters are solar disk 
source and extended emission, the likelihood is maximized over 
the other components. The analysis was performed for four dif- 
ferent estimates of the solar modulation. 




Fig. 11. EGRET Sun-centred counts maps (top to bottom)>100 
MeV, >300 MeV and 100-300 MeV. The colorbar shows the 
counts per pixel. The area is 20° on a side and the maps are 
Gaussian smoothed to 3°. 

For the different models of solar modulation, the expected 
values of inverse Compton emission are in agreement with the 



Table 4. Sum of disk and extended components of solar emission 
for different energies. Values are from the Bayesian method. 



Energy (MeV) 


Total flux (10-' crrrV 1 ) 


> 100 


(4.44±2.03) 


> 300 


(2.07±0.79) 


100 - 300 


(3.49±1.35) 



data within lcr, with rather big errors bars. This also means that, 
because of the limited sensitivity of EGRET, it is not possible to 
prove which model better describes the data. We found probabil- 
ities corresponding to 2.7, 3.6/4 and lcr above 100 MeV, above 
300 MeV (including/excluding J232 1-0328) and for 100-300 
MeV respectively for the significance of the detection. However, 
we chose the naive model of 1000 MV solar modulation in order 
to derive the spectra of the extended and disk emission. This is 
for two reasons: first this model produces the highest value of the 
likelihood; second it should be the most realistic for the period 
'91-'95, since the data cover the total period of solar maximum 
and the upcoming solar minimum. Hence, we report only the re- 
sults for the naive case of 1000 MV solar modulation. We used 
both frequentist and Bayesian approaches to analyze the results. 
The statistical methods are described in the appendix. 

Values of the best fit fluxes, lcr errors, counts and mean val- 
ues with errors are given in Table [3] Best fit values are calcu- 
lated from the log-likelihood ratio statistic, while mean values 
are from the Bayesian method. Counts are for the maximum 
likelihood values, as are the background, the 3C279 fluxes and 
the probability of the null hypothesis. Values in the first column 
were used to have more statistics for detection, while separate 
energies are best for the analysis. For 100-300 MeV range error 
bars can not be determined by the frequentist method. We also 
give the probability of the null hypothesis (i.e. zero flux from 
disk and extended emission) using the likelihood ratio statis- 
tic. The table contains two cases above 300 MeV with different 
fluxes of the background sources. The first is obtained including 
the source J232 1-0328 whose spectral index was fixed at a typ- 
ical value of 2, since it is not given in the 3EG catalogue. The 
second case does not include that source, since it has not been 
detected at energy above 300 MeV. We took the second case for 
the following analysis. These two cases give also an estimate of 
the uncertainty of our results. 

The log-likelihood ratio for the three energy ranges we ana- 
lyzed are displayed in FigQ~2]as a function of solar disk flux and 
extended flux. Colors show different values of the ratio, obtained 
by allowing the background and the QSO 3C279 component to 
vary to maximize the likelihood for each value of disk flux and 
extended component, while contour lines define 1, 2 and 3cr con- 
fidence intervals for the two separate parameters. Marginal prob- 
abilities of the two components calculated with the Bayesian 
method and cumulative probability as function of the disk and 
extended fluxes are shown in Fig. [13] The sum of the two com- 
ponents, disk and extended, is given in Table [4] Counts of the 
source components centred on the Sun resulting from the fit- 
ting technique above 300 MeV are shown in Fig[l4] A summary 
of the main results is given in Table [5] The solar disk spectral 
index we found is 2.4+yg. Errors are calculated from the uncer- 
tainties on the integrated fluxes. Figure dT~5b shows the solar disk 
spectrum obtained imposing a constant spectral index of 2.4, its 
mean value, and following a simple power law up to 2 GeV For 
the extended emission, we found a spectral index of 1.7^ g. The 
extended spectrum is shown in Fig. lfl6l l, As for the disk spec- 
trum, the spectral index has been fixed at the mean value. The 
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Table 3. Results for different energy ranges and for the naive model of 1000 MV solar modulation. Fluxes are given in cm _2 s _1 and 
intensities in crrr 2 s _1 sr _1 . For details see text. 



Model 




>100MeV 


>300 MeV with J2321-0328 


>300MeV no J2321-0328 


100-300 MeV 


OUU1LC 11 UA 


best fit 

mean 

counts 


4. (S8 +4 I / xirr 8 

t.UO 3 8? A 1U 

(5.42±3".20)xl0- 8 
28 


J.JO j g j AIL/ 

(3.03±1.10)xl0- 8 
20 


3 38 +li!8 x10- 8 

J.JO j A 1U 

(3.65±1.71)xl0- 8 
20 


1 4. 04 v 1 rr 8 

(14.2±9.1)xl0- 8 

39 


PxtpnHpH fliiY 


best fit 

mean 

counts 


3 83 +2 ' 8 x10~' 

(3.89±2.16)xl(r 7 

273 


i C4+I.U6 X 1(W 

1 . -J ~ j |g A1W 

(1.49±0.67)xl0- 7 
116 


1 7S+ 1US xlO"' 

113 

(1.70±0.87)xl0- 7 
132 


1.15 xlO - ' 

(2.07+1.34)xl0- 7 

83 


Extended solar model 


flux 


2.18 xlO"' 


0.90x10-' 


0.90 xlO-' 


1.28 xlO-' 


R a p V oToi i nH 

l J cl^lVtll WL4.11LI 


in l^Hl^i l y 

counts 


3 50x10^ 
2220 


1.13X10"" 5 
750 


1.12xl0~ 5 
741 


2.48xl0- :> 
1603 


3C279 


flux 
counts 


4.88x10-' 
103 


1.30x10-' 
29 


1.30x10-' 
29 


5.85x10-' 
79 


Prob. (flux=0) 




6.4xlO- J 


2.3xl0- 4 


9.7xl0- 5 


0.24 


Strongest bkg source 
Total bkg sources 


flux 
counts 


3.82x10-' 
113 


1.28x10-' 
27 


0.63x10-' 
21 


2.54x10-' 
66 


Moon 


flux 
counts 


5.55x10-' 
37 


0.57 xlO-' 
4 


0.57x10-' 
4 


4.98x10-' 
29 



Table 5. Fluxes used to produce the plotted solar spectra. Fluxes 
are in 10~ 7 cm^s -1 



Source 


100-300 MeV 


>300MeV 


Extended 


2.1 + 1.3 


1.7±0.9 


Model extended 


1.3 


0.9 


Disk 


1.4+0.9 


0.4±0.2 


Seckel's disk model 


0-1.1 


0.1-0.5 



spectrum is compared with the modelled spectrum for "naive" 
cases of 500 MV and 1000 MV solar modulation. 




Energy (MeV) 



Fig. 15. Solar disk spectrum. The orange regions defines the pos- 
sible values obtained by varying the mean flux within lcr errors 
and for y=2.4, the mean value of the spectral index. 



6. Tests of analysis procedure 

In order to perform a test of the solar detection, we divided 
the total exposure which we used for the analysis described 




Energy (MeV) 



Fig. 16. Solar extended spectrum. Gray regions define the pos- 
sible values obtained by varying the mean flux within lcr errors 
and for y- 1 .7, the mean value of the spectral index. Black line is 
the model for the naive case of 1000 MV modulation. For com- 
parison, the pink line shows the model for the naive case of 500 
MV modulation. 

above, into two periods with about the same exposure time. 
Then, we analyzed the data with the same fitting technique ex- 
plained above. Both these time periods contain about 3.6x1 8 
cm 2 s of exposure on the Sun. Assuming equal periods and a 
constant flux from the Sun, the two periods should give com- 
patible fit results. As an example, the analysis was performed 
above 300 MeV, where the detection was most significant and 
for the naive case and solar modulation 1000 MV. The first pe- 
riod contains only 3C279, J 1409-0745 and the moon, while in 
the second one there are only the other background sources and 
the moon; 3C279 is not present. The total emission from the Sun 
was then detected at a level of 2.6cr and 3.4c for first and second 
period respectively. The values obtained for the different emis- 
sion components are consistent within lcr. Moreover, they are 
also in agreement with those obtained for the total period. This 
confirms the validity of our method and the solar detection. 
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In order to exclude possible contamination of some faint so- 
lar flares detected by CGRO/COMPTEL during the same period, 
we performed the fitting method also with solar flares times ex- 
cluded. Values will not be reported since, as expected, we did 
not find any change in the results. 

7. Discussion 

For the extended solar emission, Fig. ( fToT i shows that our model 
is fully consistent with the measured spectrum. The main un- 
certainty in the model is the electron spectrum (the solar ra- 
diation field and the physics of the inverse Compton are pre- 
cisely known). The electron spectrum at the relevant energies 
has about a factor 2 experimental uncertainty. In fact the mea- 
sured extended fluxes are higher by a factor ~2 (although not 
very significant) which suggests our electron spectrum could be 
too low. 

Regarding the solar disk flux, the flux in all cases is in 
ag reement with the th eoretical value for pion-decay obtained 
bv lSeckel et all (119911) . The Galactic background obtained with 
the fitting technique is compatible with th e expected Galactic 
emission (GALPROP) dStrong et al.ll2004bl) . The flux of 3C 279 
of 7.2 xlQ- 1 cm "V 1 , sum of the flux for 100-300 MeV and 
above 300 MeV, is in agreement with the EGRET catalogue. 
In the catalogue for period 1 1 used for this analysis the flux is 
(7.94 ± 0.75) x lO^cm-V 1 above 100 MeV, while for period 
12 (the other viewing period used in the analysis) there is hardly 
any contribution. 

With respect to the a nalysis of EGRET data performed 
by [Thompson et al.l d!997l) . our study has some improvements 
which explain why we succeeded in detection. Instead of exclud- 
ing data near the sources 3C279 and the moon, they are included 
in our analysis. This leads to more exposure using contributions 
from all observing periods. The special dedicated analysis soft- 
ware for moving targets and the inclusion of the extended emis- 
sion from the Sun produced a more realistic prediction of the 
total data and hence a more sensitive analysis. 

8. Conclusions 

In this paper the theory of gamma-ray emission from IC scat- 
tering of solar radiation field by Galactic CR electrons has been 
given in detail. Analyzing the EGRET database, we find evi- 
dence for emission from the Sun and its vicinity. For all mod- 
els the expected values are in good agreement with the data, 
with rather big errors bars. This also means that, because of the 
limited sensitivity of EGRET, it is not possible to prove which 
model better describes the data. The spectrum of the moon has 
also been derived. 

Since the inverse Compton emission from the heliosphere is 
extended, it contributes to the whole sky foreground and has to 
be taken into account for diffuse background studies. Moreover, 
since the emission depends on the electron spectrum and its 
modulation in the heliosphere, observations in different direc- 
tions from the Sun can be used to determine the electron spec- 
trum at different position, even very close to the Sun. To distin- 
guish the modulation models (eg. with GLAST) an accuracy of 
~10% would be required. With the point source sensitivity of 
GLAST around 10~ 7 cnT 2 s _1 per day above 100 MeV, it will be 
possible to detect daily this emission, when the Sun is not close 
to the Galactic plane. In additional a precise model of the ex- 
tended emission, such as that presente d here, would be requir ed 
in searches for exotic effects such as in Fairb airn et al.l d2007l) 
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Fig. 12. Logarithm of the likelihood ratio (In (L/Lq)) as a function of the solar disk flux and an extended component for (top to 
bottom) >100 MeV, >300 MeV and 100-300 MeV . Color contours are different values of the ratio, as explained in the colorbar. 
Contours are obtained by allowing the background and the QSO 3C279 component to vary to maximize the likelihood for each 
value of disk flux and extended component. Contour lines define 1, 2 and 3cr intervals for the 2 separate parameters. 
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Fig. 13. Marginal probability of the two components calculated from eq |A.4| and cumulative probability as function of the extended 
(upper panel) and disk (lower panel) fluxes for different energy ranges, (left to right) above 100 MeV, above 300 MeV and for 
100-300 Me V 





Fig. 14. Counts of the source components resulting from the fitting technique above 300 MeV. Left to right: solar disk, inverse 
Compton, QSO 3C279, moon and background sources. Colorbars show counts per pixel. 
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Appendix A: Statistical method 

The likelihood statistic of binned data is the product of the prob- 
ability of each pixel: 

L=]~[p, 7 (A.l) 

U 

where 



Pij = JL — r ~ (A.2) 
n u \ 



is the Poisson probability of observing n, ; counts in pixel (ij) 
where the number of counts predicted by the model is 0y. For 
each case analyzed, we obtained the logarithm of the likelihood 
ratio log (L/Lq) where L is the likelihood of the data with the 
all components included, and Lq is the null hypothesis. We used 
both frequentist and Bayesian approach to analyze the results. 
Regarding the first method, -2 log (L/L m ) is distributed as x^, 
where n is the nu mber of param eters held fixed and L m is the 
global maximum (IStrongHl985h . If 2 parameters are involved, 
we used the likelihood ratio test to determine the joint signifi- 
cance detection of the components, where -2 log (L/Lq) is dis- 
tributed as x\- F° r the components separately it follows the;^ 2 
distribution. 

In order to have a better treatment of our results, we per- 
formed also a Bayesian analysis treating lik elihood as probabil - 
ity with uniform prior distribution (see eg. IStrong et al.| [2.005). 
This is independent of assumptions about distribution of the like- 
lihood ratio. Hence, the mean value of one parameter is given by: 

Xu, = J XiP(xi) dxi (A.3) 

where 



*)= f 



P(Xi)= | P(x)d"- X x (A.4) 

i*i 

is the marginal probability of x, with x^ any value with i^j, n is 
the number of the parameters and P(x) is the normalized joint 
probability distribution of x, such that 

j P(x)d"x = 1 (A.5) 
The mean square error is then given by: 

Axf = J ( Xi - x m fP{Xi)dXi (A.6) 

When we are interested in the errors in the sum of the two com- 
ponents (i.e. disk and extended emission) we have to take into 
account their covariance, since they are correlated. Hence, 

(A(xi + x k )f = (A*;) 2 + (A*,) 2 

+ 2 J(xi - x i0 )(x k - x k0 ) P(x)d n x (A.7) 



